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t is possible to extract work from a quantum-mechanical system whose dynamics 
is governed by a time-dependent cyclic Hamiltonian. An energy bath is required to 
operate such a quantum engine in place of the heat bath used to run a conventional 
classical thermodynamic heat engine. The effect of the energy bath is to maintain 
the expectation value of the system Hamiltonian during an isoenergetic process. It 
is shown that the existence of such a bath leads to equilibrium quantum states that 
maximise the von Neumann entropy. Quantum analogues of certain thermodynamic 
lO ■ relations are obtained that allow one to define the temperature of the energy bath. 



(N 



C " 3 ' Keywords: quantum-mechanical engine; von Neumann entropy; 

thermodynamic equations of state 

O " 

o 



Keywords: I 



CU 1. Introduction 



A classical thermodynamic heat engine converts heat energy into mechanical work 
by using a classical-mechanical system in which a gas expands and pushes a piston 
in a cylinder. Such a heat engine obtains its energy from a high-temperature heat 
reservoir. Some of the energy taken from this reservoir is converted to mechanical 
work. A heat engine is not perfectly efficient, so some of the energy taken from the 
heat reservoir is not converted to mechanical energy, but rather is transferred to a 
low-temperature reservoir (Planck 1927). 

A classical heat engine running between a high-temperature reservoir and a 
low-temperature reservoir achieves maximum efficiency if it is reversible. While it 
is impossible to construct a working heat engine that is perfectly reversible, in the 
early 19th century Carnot proposed a mathematical model of an ideal heat engine 
that is not only reversible but also cyclic (Carnot 1824). The Carnot engine consists 
of a cylinder of ideal gas that is alternately placed in thermal contact with high- 
temperature and low-temperature heat reservoirs whose temperatures are Tjj and 
Tc, respectively. 

Instead of a classical system of gas, we consider here a quantum-mechanical 
system of consisting a single particle in contact with a reservoir. The system is 
described by a time-dependent cyclic Hamiltonian. The statistical ensemble of many 
such identically-prepared systems that we consider here is characterised by a density 
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matrix. We assume that the system interacts weakly with its environment. It is 
known that it is possible to extract work from such a system (e.g., Geva & KoslofT 
1994; KoslofT, et al. 2000). 

In particular, if the evolution of the density matrix is reversible, then we can 
construct a quantum-mechanical Carnot engine (Bender, et al. 2000). The purpose 
of the present paper is to investigate the properties of quantum Carnot engine. 
This is of importance because it improves our understanding of the relationship 
between thermodynamics and quantum mechanics, an area of considerable interest 
in quantum theory (e.g., Leff & Rex 1990). 

This paper is organised as follows: First, we explain the concept of a quantum 
engine by generalising the two-state model considered by Bender, et al. (2000) to an 
infinite-state square-well model and derive equations of states for isoenergetic and 
adiabatic processes. These results lead naturally to the quantum analogue of the 
Clausius equality for a reversible cycle. We show that, unlike the result in classical 
thermodynamics, the Clausius relation obtained here is not based on the change of 
entropy. 

We then introduce the von Neumann entropy and obtain the maximum entropy 
state subject to isoenergetic requirements. We demonstrate that the maximum von 
Neumann entropy is consistent with the thermodynamic definition of entropy, as 
well as the equations of states for quantum Carnot cycle. As a consequence, we 
are able to determine the temperature of the energy bath directly from two of the 
diagonal components of the density matrix. We conclude by discussing several open 
problems. 



2. Quantum Carnot cycle 

We can construct a simple quantum engine using a single particle of mass m confined 
to an infinite one-dimensional square-well potential whose volume (width) is V . 
For any fixed V it is easy to solve the time independent Schrodinger equation to 
determine the energy spectrum of the system: 

E n iV) = (2.1) 

We assume that the width of the well initially is V = V\ and that the initial energy 
of the system is a fixed constant Eh- The initial state if>(x) of the system is a linear 
combination ip{x) — J2 a n4>n(x) of the energy eigenstates <j> n (x). Thus, 

oo 

J2PnEn(Vl) = E H , (2.2) 
n=l 

where p n — \a n \ 2 , and Eh is bounded below by Eh > Ei(V\). Note that the 
pure state ip{x) characterises a typical element of the ensemble, whose statistical 
property is thus determined by the density matrix, given, in the energy eigenstates, 

by Pmn — Pn^mn- 

Starting with the above initial configuration, we define a quantum cycle as fol- 
lows: First, the well expands isoenergetically; that is, the width of the well increases 
infinitely slowly while the system is kept in contact with an energy bath. Note that 
the quantum adiabatic theorem (Born & Fock 1928) states that if the system were 
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isolated during such an expansion, the system would remain in its initial state. 
That is, the absolute values of the expansion coefficients \a n \ would remain con- 
stant. Thus, if the system were isolated, the energy of the system (the expectation 
value of the Hamiltonian) would decrease like V~ 2 . 

However, during this expansion, we simultaneously pump energy into the sys- 
tem in order to compensate this decrease of the energy. Thus, the well, which can 
be viewed as a one-dimensional piston, expands in such a way that the expectation 
value of the Hamiltonian is held constant by exciting higher energy states. This 
implies that, unlike the isothermal process in classical thermodynamics, the tem- 
perature of the 'bath' in an isoenergetic process is in fact changing. The explicit 
volume dependence of temperature is derived in Section 5. 

During such an isoenergetic expansion, mechanical work is done by the force 
(one-dimensional pressure) P on the walls of the well. The contribution to this 
force from the nth energy eigenstate is 

Hence the force P is given by the expectation value P = ^2p n fn- Using the relation 
/„ = 2E n /V and Eq. (2.2), we obtain the equation of state during an isoenergetic 
process: 

PV = 2E H . (2.4) 

This result is identical to the corresponding equation of state for an isothermal 
process of a classical ideal gas, if we make the identification 2E H <-> kT H . Note 
that the expansion coefficients a n of the wave function change as a function of the 
width V while the well expands isoenergetically from V\ to V 2 - 

Second, the system expands adiabatically. During an adiabatic process the eigen- 
states 4> n {x) change as a function of V, but the values \a n \ remain constant. There- 
fore, the expectation value of the Hamiltonian E = J2p n E n (V) decreases during 
the process because each E n decreases with increasing V while all p n are kept fixed. 
The force P in this case is determined by differentiating the energy E with respect 
to V. Thus, the equation of state during an adiabatic process for the square- well 
potential is 

PV 3 = 2V 2 E H , (2.5) 

which is a quantum analogue of the corresponding equation for a classical ideal 
gas. The system expands adiabatically until its volume reaches V — V3. At this 
point the expectation of the Hamiltonian decreases to Ec- Because the squared 
coefficients p n of the wave function remain constant during an adiabatic process, 
the value of E c is given by E c = {V% /V 2 )E H . 

Following the adiabatic expansion, the system is compressed isoenergetically 
until V — V4, with the expectation value of the Hamiltonian fixed at Ec, and 
finally it is compressed adiabatically until the width of the system returns to its 
initial value V — V\. This cycle is reversible, and we find that the efficiency of the 
quantum engine is given by 

f = l-|£. ( 2 -6) 
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a formula analogous to the classical thermodynamic result 77 = 1 — Tq /Th of Carnot 
(1824). 

3. Quantum Clausius relations 

We have demonstrated above an example of a quantum system with a slowly- 
changing time-dependent Hamiltonian from which work can be extracted. The key 
concept introduced here is an energy bath that maintains the expectation value 
of the Hamiltonian. This idea can be developed further to establish the quantum 
counterparts of the classical thermodynamic relations. To begin with, let us first 
consider the theorem of Clausius. 

During an isoenergetic expansion, the amount of energy transferred to the sys- 
tem to maintain the expectation of the Hamiltonian is determined by the integral 

rVl 

Qh = dVP(V) = 2E H ln(V 2 /V 1 ), (3.1) 

where P(V) — 2EhV~ 1 . The amount of energy absorbed during an isoenergetic 
compression can be determined similarly with the result Qc — — 2Eq ln(Va/V4). 
Thus, for a closed, reversible cycle we obtain the quantum Clausius equality 

Q± + 9<L = 0, (3.2) 

because V2/V1 = V3/V4 for a closed cycle. 

In classical thermodynamics the Clausius equality states that the total change 
§ dS of entropy in a reversible cycle is zero, where the differential of entropy 
dS = dQ/T is the ratio of the absorbed heat and the bath temperature. There- 
fore, the relation (3.2) suggests that in quantum mechanics the entropy change in 
an isoenergetic process might be given by dQ/E, the ratio of the absorbed energy 
to the bath energy. However, as we show below, this quantity does not determine 
the amount of entropy change, and hence the quantum Clausius relation is not a 
condition for entropy. Instead, it merely implies that the total energy absorbed, for 
given bath energies, must vanish in a reversible cycle. If the cycle is irreversible, 
then we have § dQ/E < 0. 

4. Entropy maximisation 

To understand the change of entropy in a quantum Carnot cycle we consider the 
von Neumann entropy 

00 

S = -^p n \np n (4.1) 

n=l 

associated with the density matrix p mni expressed in terms of the energy eigen- 
states. Let us first consider the isoenergetic process and introduce a scaling param- 
eter A such that V = XVi . Then, during an isoenergetic process the probability p n 
that the system is in the nth eigenstate changes in A, subject to the constraints 

00 00 

^p n = l and ^n 2 p„ = A 2 , (4.2) 

n—1 n—1 
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where we have chosen units in which Tr 2 h 2 /(2m) = 1 and taken the initial condi- 
tion Eh — Ei(Vi). There are infinitely many combinations of p n satisfying these 
constraints, each of which determines the entropy. Hence, quantum mechanically 
it appears that any one of such states is allowed in an isoenergetic process. How- 
ever, there is a unique density matrix that satisfies the thermodynamic requirement 
of the change of entropy, and this is the state that maximises the von Neumann 
entropy. 

To show this, we must determine the density matrix that maximises the entropy 
subject to the constraints (4.2) and then obtain the associated entropy S. We can 
perform this maximisation directly without using Lagrange multipliers. The two 
constraints in (4.2) imply that two of the diagonal components, say pk and pi, of 
the density matrix are determined. Therefore, if we differentiate the constraints 
(4.2) with respect to p n , we find that 

dpk dpi 2 dp k 2 dpi 2 

|-__ = -l and k z - hi 7- = -n • 4.3) 

Op n Op n dp n 0p n 

Solving these linear equations, we obtain 

dpk _ n 2 — I 2 dpi _ n 2 - k 2 

0p n ~ k 2 ~l 2 d Pn ~ l 2 -k 2 [ ] 

for all k,l^n. Substituting these expressions into dS/dp n = and solving for p n , 
we get 

which is also valid for any choice oik ^l. 

Given the maximum entropy state (4.5), we can now reexpress the constraints 
in the form 

A 2 f>" 2 =f>v* 2 , ^ 

n— 1 n— 1 

where we have defined a = (pk/pi) 1 ^^ ~' 2 '- This is obtained by substituting (4.5) 
in (4.2), solving for pi, and equating the resulting expressions. Similarly, using the 
parameter a, the state represented in (4.5) can be expressed as 

a" 2 

Therefore, the corresponding von Neumann entropy is 

S = {I 2 - X 2 )\na-\npi. (4.8) 
This is independent of I because of the identity 

p k a~ h2 =pia~ l \ (4.9) 
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Thus, given the width scale A = V/V\ of the potential well, the relation (4.6) 
determines the value of a, from which we can determine both p n and S as functions 
of the single length scale parameter A. 

We can now analyse the entropy associated with the isoenergetic process. We 
note first that the change of entropy associated with (4.8) does not agree with the 
quantity dQ/E. More specifically, we verify that only asymptotically for large values 
of A, that is, large energies, the two quantities are proportional to each other. In 
particular, because the value of S in (4.8) is independent of the choices of k and I, 
we can set k = 2 and I = 1 without loss of generality. Then, asymptotically we have 
Pi ~ P2 for large A. Furthermore, the maximum entropy state {p n } satisfies pk > pi 
for all k < I. Therefore, a ~ 1 and for large values of A we use the asymptotic 
relation 

n=l v 

to determine the entropy change. The result gives AS ~ In A when the width 
changes from V\ to V = XVi. On the other hand, from (3.1) we deduce that 
dQ/E — 2 In A for any value of A > 1. This establishes the asymptotic propor- 
tionality. However, dQ/E does not in general determine the entropy change. 



5. Equilibrium-state temperature 

Next, we show that the maximum von Neumann entropy (4.8) is consistent with 
thermodynamic considerations. Recall that an isoenergetic process must proceed 
infinitely slowly in order that the energy bath at any point during the process be 
in thermal equilibrium. Thus, this process requires an infinite amount of time to 
be realised. Then from the principles of statistical mechanics (Schrodinger 1952) 
we deduce that, if we set V\ = 1, the temperature of the bath is given by T = 
— l/(A 2 lna). This is because we can make the identification a n — cnp(—E n /T) 
in Eq. (4.7). On the other hand, temperature is given in thermodynamics by the 
formula 

i-^ (51) 

where dQ = PdV. In the present consideration we have dQ — 2\~ 1 d\ from the 
equation of state PV = 2Eh and the relation dQ = PdV, while from (4.8) we find 
that dS/dX = — 2Alna. Therefore, the thermodynamic definition of temperature in 
(5.1) gives 

T = -T2T- > ( 5 ' 2 ) 
X z ma 

which agrees with the statistical mechanical consideration above. We conclude 
therefore that the maximum von Neumann entropy determines the density ma- 
trix of the system in an isoenergetic process. Furthermore, the temperature of the 
energy bath can be obtained from (5.2), where a is determined by any two diagonal 
components pk, pi of the density matrix. In particular, the temperature of the bath 
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Figure 1. 


Phase diagram in the S 


— A plane. The entropy of the system increases during an 



isoenergetic expansion, satisfying the constraint {E n ) = Eh, then remains constant for an 
adiabatic expansion, decreases for isoenergetic compression with the constraint {E n } — Ec 
to its initial value, and remains constant for adiabatic compression. 

changes continuously in A during an isoenergetic process, even though the energy 
of the system is held fixed. 

We now consider the adiabatic process. In this case, there is no net energy 
absorbed by the system, so the thermodynamic entropy remains constant. On the 
other hand, because the diagonal components p n of the density matrix are constant 
during an adiabatic process, the von Neumann entropy (4.1) also remains constant, 
in agreement with the thermodynamic consideration. 

These results can be expressed by the phase diagram in the entropy-volume 
plane, as shown in figure 1. The diagram shows that the thermodynamic constraint 
V2/V1 — V3/V4 is satisfied by, and only by, the maximum von Neumann entropy. 
Thus, we emphasise that the quantum state corresponding to the maximum von 
Neumann entropy, as obtained by microscopic analysis, is entirely consistent with 
the thermodynamic equations of states. 

6. Discussion 

In the analysis above we have considered only the square-well potential model, 
which can be interpreted as the quantum analogue of the classical ideal gas. A 
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natural extension of this work would be to consider other Hamiltonians. For ex- 
ample, for a harmonic potential whose energy eigenvalues are E n = (n + ^)hu), a 
dimensional argument shows that the characteristic length scale is given by 1/y/uj. 
Therefore, the equations of states for the harmonic potential become identical to 
those for the square- well potential, although the volume dependence of the entropy 
is different from (4.8). For other Hamiltonians, however, the equations of states are 
not in general identical to those obtained here, and it would be interesting to find 
explicit examples of other models exhibiting nonideal behaviours. 

Another important issue is to understand whether there is any role played by 
the geometric phases, a concept that does not have an analogue in the classical 
Carnot cycle. Although the quantum Carnot engine is indeed cyclic in terms of p n , 
the coefficients a n of the wave function of any one of the representatives in the 
ensemble inevitably pick up geometric phases as the system goes through a cycle. 
Therefore, strictly speaking, the wave function is not cyclic in the quantum Carnot 
cycle. The question then is whether there is any physically observable evidence of 
the geometric phase in the present context, and if so, how do we represent that in 
the mixed state context (cf. Uhlmann 1996). 

Finally, another interesting idea that should be studied is the possibility of con- 
structing a cyclic quantum engine (or a quantum refrigerator by reversing the cycle) 
that requires only a finite amount of time to complete. Note that each cycle of a clas- 
sical Carnot engine is infinitely long. However, in quantum mechanics, it is known, 
for example, that by choosing a specific form of time dependent Hamiltonians, it is 
possible to construct a quantum-mechanical cycle in an essentially arbitrary short 
time scale (Miclnik 1986). It would be of interest to determine if such an idea can 
be applied to create a finite-time Carnot engine, or whether such dynamics would 
be incompatible with the requirement of maximising the von Neumann entropy. An 
analysis on the performance of finite-time cycle in the presence of friction has been 
studied in Feldmann & Kosloff (2000). 
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